Coarse Grained Microscopic Model of Glass Formers 
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We introduce a coarse grained model for atomic glass formers. Its elements are physically mo- 
tivated local microscopic dynamical rules parameterized by observables. Results of the model are 
established and used to interpret the measured behaviors of supercooled fluids approaching glass 
transitions. The model predicts the presence of a crossover from hierarchical super- Arrhenius dy- 
namics at short length scales to diffusive Arrhenius dynamics at large length scales. This prediction 
distinguishes our model from other theories of glass formers, and can to be tested by experiment. 



Introduction. This article introduces a microscopic 
model of glass forming supercooled liquids. It is a coarse 
grained model with local dynamical rules, the elements 
of which can be measured. The rules contain only the 
physical features of detailed balance, dynamic facilita- 
tion, and persistence of particle flow direction. Our anal- 
ysis of the model indicates that these features are all 
that is required to explain the time dependent and ther- 
mal behaviors of supercooled glass forming liquids. For 
such systems, relaxation times grow rapidly with decreas- 
ing temperature T, usually more rapidly than exponen- 
tially in T _1 . This non- Arrhenius (or super- Arrhenius) 
temperature dependence often coincides with a relatively 
large change in heat capacity at the glass transition. On 
the other hand, liquids whose relaxation times grow in 
an Arrhenius fashion with decreasing temperature often 
display small changes in heat capacity at the glass tran- 
sition. Apparent correlations like these between thermo- 
dynamic properties and relaxation times have led to a 
commonly held opinion that a finite temperature thermo- 
dynamic transition somehow underlies a glass transition 
(for a review of this and other views, see fl]). In con- 
trast, our model explains the empirically observed cor- 
relations of relaxation and calorimetric properties at the 
glass transition without recourse to unusual or precipi- 
tous thermodynamics. 

Our model is described in terms of a grid or lattice 
in space-time. Underlying this discrete picture are cer- 
tain ideas about continuous atomic motions in a super- 
cooled liquid. The first idea is that particle mobility is 
sparse. Most atomic motions are small amplitude vibra- 
tions and not diffusion steps. Most molecules are effec- 
tively jammed. The second idea is that particle mobility 
in a glass former is the result of facilitation H0. Specif- 
ically, a region of jammed atoms can become unjammcd 
and thus exhibit mobility only when it is adjacent to a 
region that is already unjammed. This idea implies that 
mobile or unjammed particles form chains in space-time. 
Chains of mobile particles have been revealed byatom- 
istic computer simulations of supercooled liquids [4[ • The 
third idea is that mobility carries a direction, and it is 
possible for this direction to persist in a glass former for 
significant periods of time. This convective-like quality 
of particle motion is also revealed by the atomistic com- 
puter simulations of supercooled liquids Q . 



To help illustrate how these ideas motivate a lattice 
model, Fig. 1(a) depicts a mostly jammed fluid of spheres 
(drawn in two dimensions for artistic convenience). It 
shows, schematically, how mobility is found in unjammed 
regions, and that directed motion in one of those regions 
might make similarly directed motion possible in adja- 
cent regions. The arrows are drawn pointing parallel to 
the mean direction of facilitation, and anti-parallel to the 
mean direction of particle motion. Given that particles 
singly occupy finite volumes, the possible influence of fa- 
cilitation falls in a cone, as indicated by the dashed lines 
in Fig. 1(a). 

The point of each cone in Fig. 1(a) depicts a position 
that is made accessible to an atom by prior motion in the 
vicinity of that point. As large voids dissipate quickly in 
a dense material, the most likely motions are collective 
and complicated shuffling of several particles. This shuf- 
fling could have a number of causes. One possibility is 
the presence of local density that is slightly less than 
a constraining value. Whatever its precise nature, the 
collective motion will have a definite direction for a pe- 
riod of time. In highly networked systems, directional 
persistence will quickly dissipate, decreasing its effect on 
subsequent motion. 

To capture these physical features, we use the coarse 
grained description depicted in Fig. 1(b). Space is di- 
vided into square or cubic regions with side lengths com- 
parable to, but not smaller than, the bulk correlation 
length of the fluid. As such, equilibrium fluctuations in 
one cell are uncorrelated with those in another at the 
same point in time. Further, time is divided into dis- 
crete intervals or time steps of length St long enough to 
discriminate whether mobility (i.e., significant particle 
displacement) has or has not occurred in a particular cell 
during that interval. An empty cell in Fig. 1 (b) indicates 
a region that exhibits no mobility during that time step. 
A cell containing an arrow indicates a region that does 
exhibit mobility during that time step, and the direction 
of the arrow is anti-parallel to the direction of motion 
and therefore indicates the direction of facilitation. 

This coarse grained catalog or description of dynami- 
cal behavior can be constructed from the trajectories of 
any system. In this sense, there is no approximation in- 
volved in its use. Were the system an ordinary liquid, 
a typical equilibrium configuration of the coarse grained 
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FIG. 1: (a) Facilitation of motion in a nearly arrested fluid. 
The arrows indicate direction of facilitation, anti parallel to 
direction of flow. Motion is facilitated in cones indicated by 
dotted lines. The white circles depict positions of atoms at 
one point in time; each grey circle, near a white circle of 
the same size, indicates the position of that near atom at a 
later point in time, (b) Coarse grained lattice description. 
The excitation at the bottom facilitates cells a and b, the 
excitation at the top, cells c and d. (c) Possible series of 
steps leading to a diffusive step in d = 2. 



model would have essentially all regions occupied. For 
a sufficiently supercooled liquid, however, the concentra- 
tion of arrows will be low. In that case, the thermody- 
namics of the material should be consistent with that of 
non-interacting arrows on a lattice. As a simplification, 
we shall consider only discrete orientations of the arrows, 
restricting them to lie along any one of the principal di- 
agonals of its cell. Placing arrows along cell diagonals 
rather than edges allows for a cone of influence with pos- 
itive angle. In particular, an arrow at a specific cell and 
time step can facilitate the creation or destruction of an 
arrow in a nearest neighbor cell at the next time step, 
where the vector drawn from the specific cell towards 
that nearest neighbor has a positive inner product with 
the arrow. For example, in Fig. 1(b), the arrow at the 
lower left can facilitate the birth of an arrow at the next 
time step in either cell labeled a or b. In this case, the 
cone of influence spans a quadrant in the square lattice, 
or an octant in the cubic lattice. Different cones of in- 
fluence can be defined with different choices of spatial 
grids. 

The model. As the preceding discussion indicates, our 
model describes the dynamics of the vector field n(x) = 
n(x)v(x), where n(x) — 1 or 0, depending on whether 
atoms in cell x are mobile or immobile, respectively, and 
v(x) is the unit vector pointing in the direction of the ar- 
row in cell x. For a square (d = 2) lattice, v(x) will be one 
of the four vectors (±1,±1) /\/2. The equilibrium con- 
centration of mobile cells is denoted by c — (n(x)) (taking 
the lattice spacing to be unit length). The pointed brack- 
ets indicate equilibrium ensemble average. The system is 
isotropic, so that (v(a ? )) = 0, and only trivial equal-time 
correlations exist so that (n(x)n(y)) = for x ^ y . Ac- 
cordingly, the equilibrium distribution for the vector field 



is P({n(x)}) =T[ S p[n(x)}, where 
p{n)=g- n {l-ct 



(1) 



Here, g is the number of equally likely equilibrium orien- 
tations of an arrow [i.e., a vector v(af)]. For the square 
(cubic) lattice, g = 4 (g = 8). 

The passage from one time interval, t, to the next, 
t + St, is governed in our model by the following single 
site Monte Carlo dynamics. Consider the transition at 
site x of the vector field n(x) = n(x)v(x) from an un- 
excited state, n(x) = 0, to an excited state n(x) = 1, 
v(x) = w, and the corresponding backward transition. 
The dynamic facilitation described above is implemented 
by defining the corresponding rates as: 



C a [w] cjg 



n(x) = 



n(x) = w, 



(2) 



C 3 [w] (1-c) 



where the kinetic constraint is implemented by the indi- 
cated state dependent rates, with Cg[w] given by 



C*[w] = / 
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Here, indicates the i-th Cartesian component of w, 
nearest neighbors are denoted by (y, x). Since [w] docs 
not depend on the state of the cell making the transition, 
but only on the state of its neighbours, it disappears from 
the ratio between forward and backward rates, and the 
dynamics defined by (J2J satisfies detailed balance with 
respect to the distribution Q. The parameter / e [0, 1] 
in Eq. © determines the probability p(f) that a newly 
created arrow will lie parallel to its facilitating arrow, 
p(f) = [1 + f{g — l)] -1 - One of the extremes, / = 1, 
corresponds to the case where an excitation pointing in 
any direction can be created or destroyed as long as its 
site is facilitated by a neighbour. In the other extreme, 
/ = 0, only excitations pointing in the same direction as 
the facilitating neighbour can be created or destroyed. 

The above dynamics satisfies time reversal symmetry 
and detailed balance. It is designed to be physically real- 
istic for c small, the realm of dynamical arrest in a glass 
former. When c is not small, the dynamics of a realistic 
system will proceed through other mechanisms in addi- 
tion to facilitation. In general, even for c small, space 
and time coarse graining of an atomistic model will lead 
to transition rules with non-Markovian elements. We as- 
sume these elements are unimportant. 

The parameters c and / have clear physical meaning. 
Observations of trajectories can determine the probabil- 
ity, c, that a microscopic region exhibits mobility over 
a specified microscopic time period. They can further 
determine / from the directional persistence probability 
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over a similar time period. For systems where directional- 
ity in particle motion can persist over a microscopic time 
St (e.g., a nearly jammed fluid of hard spheres), / will be 
very small. For systems where directionality in particle 
motion will dissipate over that time (e.g., a highly net- 
worked fluid), / will be close to 1. One may use atomistic 
computer simulations of supercooled fluids to make the 
observations that determine c and /. Alternatively, one 
may perform microscopy on colloidal glass formers. In 
these cases, observed time scales relative to those of an 
ordinary liquid are modest in length. For glass forming 
supercooled liquids exhibiting long relative times scales, 
experiments using single molecule spectroscopic probes 
could be used. 

The parameters c and / have energetic and entropic 
contributions. Since excitations in different sites are 
statically uncorrelated, their average concentration c 
must have the usual Boltzmann temperature dependence. 
Therefore, to the extent that c is small and mainly en- 
ergetic, — In (c/g) is proportional to reciprocal tempera- 
ture, 1/T. Similarly, since / is a local probability inde- 
pendent of the state of the system, it is natural to write 
/ = ac b where a pertains mainly to its entropic contribu- 
tion, and b pertains mainly to its energetic contribution. 

Properties of the model. It is useful to begin by con- 
sidering the two extreme limits of the model, the / = 1 
case and the / = case. 

• / = 1. Here, motion of an arrow can proceed dif- 
fusively as illustrated with the seven step sequence con- 
figurations shown in Fig. 1(c). This sequence, involv- 
ing single activation on three occasions, produces a sin- 
gle diffusive step in either d = 2 or 3. As such, the 
relaxation time for a length scale I, r/ = i(i), will obey 
Tf = i(l) ~ 8tD~ l P 1 , where the dynamical exponent, Z\, 
will be a constant, while the diffusion constant, D, will be 
proportional to c for d = 2 and 3. (The quantity D is the 
diffusion constant for excitations, which is distinct from 
the diffusion constant for particle motion.) Since the dis- 
tance between arrows at equilibrium is c" 1 ^, the system 
approaches equilibrium with time scale r/=i ~ 5tc~ Al 
where for d = 2 and 3, Ai = 1 + z\/d. From numerical 
simulations of the model on a square (d — 2) and cubic 
(d = 3) lattices, we find that Ai w 2. Simulation results 
for the relaxation of the model from a disordered state 
on a square lattice are shown in the left panel of Fig. 
2. As expected from the diffusive nature of its motion, 
the behavior of this model is like that of the FA (i.e., 
one-spin- facilitated Fredrickson- Andersen 0, |(|) model. 
The Arrhenius dependence of Tf—\ on c is often called 
"strong" 0. 

• / = 0. Here, the motion is hierarchical 0,13, like 
that of the East model [aOl or its generalization to higher 
dimension [T^j, and it is crgodic for d > 2. An arrow 
can relax only through pathways involving like pointing 
arrows. For d = 2 the system is composed of a mixture 
of arrows of four different directions, with each direction 
uncoupled from the others. The dynamical behavior of 



each of the four independent components is like that of a 
North-or-East model |T(| which in turn behaves like the 
East model. In this case, therefore, t/ = o(1) ~ St l z ° where 
now the dynamical exponent is not constant but rather 
a function of c, namely zq ~ ln(<?/c)/ In 2, where c/g is 
the concentration of like directed arrows. Hence, close to 
equilibrium, the overall time scale for relaxation is t/ = o ~ 
St c~ A °, where Ao = Zo/d. This scaling is consistent with 
the simulation results for d = 2 graphed in the right panel 
of Fig. 2. The super-Arrhenius dependence on c is often 
called "fragile" Q. 

• < / < 1. Here, we adopt an interpolation formula 
by adding rates in series: t _1 (I) « (1 — /) tJ/z q (I) + 
f 7y_j (I). As such, for small /, 

t(1) « Stl Zo (l + D f . (4) 

The corresponding equilibrium relaxation time is there- 
fore 

r«<5ic- A ° + / C -( A °- Al >)~\ (5) 

These equations exhibit a crossover behavior. In par- 
ticular, for a given c and /, there is a crossover length, 
lit, such that for I > £*, relaxation is diffusive and thus 
Arrhenius, and for I < Z„, relaxation is hierarchical and 
super-Arrhenius. Thus, at large enough length scales, 
the relaxation of all systems, no matter how fragile, will 
obey diffusive behavior. The existence of the crossover 
length, I*, could be demonstrated experimentally from 
the single molecule spectroscopy of rotating nanoparticle 
probes of various lengths. The existence of a crossover 
length also implies that there is a crossover c* (or if 
one associates c with T as suggested earlier) such that for 
c < (or T < T*), the equilibrium relaxation time of 
the model is Arrhenius, and for c > c* (or T > T*), it is 
super-Arrhenius. Implicit in these remarks is an appre- 
ciation for the presence of dynamic heterogeneity |g, 1 1 1| , 
in this case including flow direction, but for limitations 
of space, we cannot discuss it further here. 

Crossover behavior was already demonstrated for what 
can be considered a mean field version of our model |12| . 
The presence of a crossover has also been argued for on 
experimental grounds ^| ■ Figure 3 illustrates the conse- 
quent temperature dependence for the possible crossover 
behaviors predicted from our model. Graphs on the left 
are drawn as Arrhenius plots. Graphs on the right show 
the same quantities but scaled in the fashion proposed 
by Angell [l| . Notice the possibility that relaxation times 
may fall above the Arrhenius line in the Angell plot. This 
specific behavior has not been anticipated. It arises in 
systems where / is both small and mainly entropic (i.e., 
independent of c) . 

Comparison with experiment. The fundamental 
control parameter is c, and while this quantity could be 
measured directly, most existing experimental data ob- 
serves the effects of changing c by changing temperature 
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FIG. 2: Relaxation of the concentration c(t) after a quench 
from an initial random state to final equilibrium concentration 
c in the d — 2 version of the model. Left panel: / = 1. Right 
panel: f — 0. Simulations were performed using continuous 
time Monte-Carlo ;5|, on system sizes L = 10 5 and averaged 
over 50 samples. Error bars are comparable to line thickness. 



at the fixed pressure of 1 atm. To compare with that 
data, we therefore use 

where cr is the concentration of mobile regions at a ref- 
erence temperature Tr, and J is the enthalpy difference 
between a mobile and immobile region. Equation JSJ is 
again based on the idea that since excitations at differ- 
ent spatial points (at equal times) are statically uncor- 
rected their average concentration c must have Boltz- 
mann temperature dependence. Since mobile regions 
should be liquid like, and immobile regions should be 
solid like, we expect that J/T should be of the order of a 
few entropy units times the typical number of molecules 
in a cell (~ 10). With lnc linear in 1/T, lnr/ = o is 
a quadratic polynomial in T^ 1 14], while lnr/ = i is a 
linear one. A useful reference point, which allows for 
a direct comparison with previous analysis (eg., Ref. 

is Tr = T1/2, the temperature at which the log- 
arithm of the relaxation time is one half that at the 
glass transition. In this case tr = t(cr) is of the or- 
der of microseconds (/xs). A common measure of fragility 
is i*i/2 = 2Tg/T]/o — 1, where T g is the glass transi- 
tion temperature |15| . A second common measure is the 
steepness index, to = d\ogT(T g )/d(T g /T) 0, i.e., the 
slope at T g in the Angell plot. If the relaxation from 
Tx/2 to T g follows the super- Arrhenius form of t/ = o, 
then F1/2 and the elementary timescale St determine 
TO: m/=0 = 32(1 - F 1/2 )- x (l + Y)(l + 2Y)~\ where 

Y = \ag((j&/6t) [1 + V 1 + 4/log(/xs/(ft)]/8. In the case 
of Arrhenius behaviour, or of an energetic crossover that 
takes place at T > T g , the steepness index is given by 
m/=i = 16(1 — F 1 / 2 )~ 1 , independent of St. 

Figure 4 shows that the ranges of fragility measures 
[TEl llq are well understood with our model. The dashed 
line corresponds to the value of m/=o in the case where 




ln(c /c) ln(c /c) / ln(c /c ) 



FIG. 3: Left panel: Arrhenius-like plot of relaxation time 
against equilibrium concentration, for various values of / = 
ac b . The cases a — 1, b = and a = correspond to purely 
Arrhenius and super- Arrhenius behaviour, respectively. The 
case = 6 = 1 leads to an energetic crossover, while a = 
.01,6 = 0, to a mainly entropic one. Right panel: same as 
before now presented as an Angell-like plot, with the glass 
transition chosen at r ~ 10 10 <5i. c x = 8/9 is the random 
concentration value for the cubic lattice. 



St ~ (is, m/ =0 = 32(1 - F 1/2 )- 1 . The full line is 
mf—x- The region between this two lines, region I in 
Fig. 4, corresponds to super- Arrhenius behaviour Tt—o 
with St < fj,s. Most glass forming liquids fall in this 
region. Liquids on the full line are either those with 
purely Arrhenius behaviour or those that undergo an 
energetic crossover. Region II below the full line, corre- 
sponds to liquids which experience the entropic crossover 
noted above. Salol (SL) and a-phenyl-o-cresol (APOC) 
fall in this region. Region III above the dashed line cor- 
responds to values of m which are too large to account 
with the present approach (except with St ^> fis, which 
seems unphysical). Since selenium forms polymer chains, 
a feature not accounted for in our model, it is perhaps 
not surprising that it falls in this region |17|. 

In Fig. 5 we show Angell plots in the regime near 
T g for five representative liquids from Fig. 4 (data from 
It shows fits to the dielectric relaxation timescale of 
3-bromopentane (3BP) and viscosity of ortho-terphcnyl 
(OTP), using tj = q, and an Arrhenius fit to the viscosity 
data of GeC>2, using r/ = i. It also illustrates the crossover 
in the the viscosity data of APOC and SL. The fits to 
the kinetic data determine the parameters of the model. 
From the coefficients of the terms quadratic in T g /T in 
logr/ = 0i an d from the linear one in logr/ = i, we obtain 
the values of J/T g = (16.7, 26.7, 29.3, 31.8, 22.6) for (3BP, 
OTP, APOC, SL, Ge0 2 ). These values of J/T g obtained 
from fitting kinetic data agree with the expectations of 
our model. 

Further, from the linear and constant terms in log t/ = o, 
we determine log(St//j,s) = (-3.3,-4.4,-4.7,-4.8) and 
log(c g / 5 ) = (-3.2,-3.4,-3.6,-3.9), for (3BP, OTP, 
APOC, SL), respectively, so that in these cases St is 
of the order of fractions of ns. These are physically 
reasonable times to discern whether or not a micro- 
scopic region exhibits mobility. We also obtain a re- 
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FIG. 4: Steepness index m against kinetic fragility Fi/ 2 for 
various supercooled liquids. 




FIG. 5: Angell plots of dielectric relaxation and viscosity. In 
the left panel, solid lines for 3BP and OTP refer to r/ = o, and 
solid line for GeC>2 refers to Tf = \. In right panel full line 
corresponds to Eq. 10 ; dotted and dashed lines indicate 17=0 
and 77=1 behaviour before and after the crossover. Symbols 
indicate experimental data. 



lation between the concentration at T\/ 2 , to that at 
the glass transition, log{c g (l — c 1 / 2 )/[c 1 / 2 (l — c g )]} = 
(-1.5,-1.4,-1.5,-1.8,-4.6), for (3BP, OTP, APOC, 
SL,Ge0 2 ). This list of ratios could in principle be 
checked experimentally. Fits to kinetic data do not de- 
termine the value of the parameter g. If we assume g « 8 
as in the lattice model, then we obtain for the previous 
liquids: logc g = (—2.3, —2.5, —2.7, —3.0, —3.6), where in 
the case of Ge0 2 we have taken ci/ 2 ~ 0(1). Notice that 
the mobility concentration at T g of the two purely super- 
Arrhenius liquids, 3BP and OTP, is at least an order 
of magnitude larger than that of the purely Arrhenius 
one, Ge0 2 . Moreover, the two liquids which display a 
crossover have intermediate values of c g . 

The order of magnitude difference between fragile and 
strong liquids (OTP and 3BP versus Ge0 2 in our ex- 
ample) in concentration of mobility excitations at T g is 
a crucial theoretical prediction of our approach. In the 
context of the model this difference has a clear origin. 
Relaxation time grows with decreasing c much faster in 
the fragile case (/ = 0) than it in the strong case (/ = 1). 
Thus, an arbitrarily specified long relaxation time (e.g., 
some minutes) coinciding to a glass transition is reached 
at a value of c in the fragile case that is much larger than 
that in the strong case. 

The phenomenological consequences of this difference 
in concentrations are important. For example, l/c 1 / d (T) 
sets the dynamic correlation length, which in turn gives 
an upper bound for the size of dynamical heterogeneities 
Assuming that each cell contains about 5 to 10 
molecules, so that Sx « 1-1.5 nm, we estimate that the 
size of dynamic heterogeneities for 3BP, OTP and Ge0 2 
are bounded from above by 5-8 nm, 6-9 nm and 15-24 
nm, respectively. As we have stressed 8], dynamic het- 
erogeneity is self-similar over a range of lengths. There- 
fore, at a specific thermodynamic state point, no single 



length scale can be attributed to dynamic heterogeneity. 
Nevertheless, experimentalists have reported results for a 
few frag ile liquids as if there is but a single specific length 
[Til |20| . The reported experimental lengths are consis- 
tent with our predicted bounds, though a definite com- 
parison between theory and experiment must await more 
complete experimental information on the full range of 
heterogeneity lengths. Since l/c g is larger in the strong 
case than the fragile case, our model predicts that the 
size of dynamical heterogeneous regions in strong liquids 
should be appreciably larger than those in fragile liquids. 
This predicted trend is opposite of that made by other 
theoretical approaches, like frustrated limited domains 
|2lj| or random first order transitions [Tflj . 

From the concentrations c g and the corresponding val- 
ues of J/Xg we are able to predict the jumps in specific 
heats at T g . According to our model, this heat capacity 
(per mol) relative to that of the solid is 

Ac p (T g )*(J/T e fc g M+0(cl), (7) 

where Af is the number of molecules that contribute 
to enthalpy fluctuations per mobile cell. With the val- 
ues obtained from the fits to the kinetic data for OTP, 
3BP and Ge0 2 , for example, we obtain, Ac p {T g )/N ~ 
2.3, 1.3,0.13, respectively. These values are in the ratio 
17 : 10 : 1, which compares well to the experimental 
values Ac p {T g )/k B ~ 13.6,9.2,0.84 22], whose ratio is 
16 : 11 : 1. 

Further comparisons between this theory and exper- 
iment illustrating subtle counter examples to the usual 
correlation between super- Arrhenius behaviour and Ac p 
[23| are left to future papers. Along with comparisons 
to equilibrium properties, illustrated herein, our model 
can be applied to out-of-equilibrium behaviors without 
further assumptions. 
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